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Abstract 

We present the results of fitting elastic pp differential cross section data 
at 23.5 < < 62.5 GeV with a novel analytic parametrization for the 
scattering amplitude. Making use of a fitting method, the errors from the 
free parameters are propagated to the imaginary part of the eikonal in the 
momentum transfer space. A novel systematic study of the effects coming 
from data at large momentum transfer is also performed. We find statistical 
evidence for the existence of eikonal zeros in the interval of momentum transfer 
5-9 GeV 2 . 
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Several authors have investigated elastic pp scattering at high energies in the context of 
the impact parameter formalism and/or eikonal approximation. Model-independent anal- 
yses, performed through fits to differential cross section data, play an important role as a 
source of empirical information for theoretical developments. 

Inspired by some of these mo del- independent approaches and based on experimental 
information presently available J7||| we retake two crucial aspects of the empirical analyses: 

(a) a systematic study of the effect of experimental informations at large momentum transfer; 

(b) development of a statistical procedure in order to estimate the eikonal uncertainties in the 
momentum transfer space. In this work we briefly report some results that can contribute 
with a better understanding of the above aspects, and bring important insights for the 
development of theoretical approaches. 

We shall analyze experimental data on elastic pp differential cross section and the p 
parameter (ratio of the forward real to imaginary part of the scattering amplitude) between 
\fs = 23.5 GeV and = 62.5 GeV. Our basic set corresponds to the data compiled and 
normalized by Amaldi and Schubert at the energies \/s = 23.5, 30.7, 44.7, 52.8, and 62.5 
GeV and momentum transfer in the interval 0.01 GeV 2 < q 2 < 9.8 GeV 2 |7|||. However, 
based on the evidence that these differential cross section data do not depend on the energy 
for momentum transfer above q 2 ~ 3.0 GeV 2 P,[l0| we shall also take account of the data 



available at ^/s = 27.5 GeV, in the interval 5.5 < q 2 < 14.2 GeV 2 0. 

Our strategy in performing a systematic study of this experimental information is to 
consider two different ensembles of data, characterized and denoted by the following. 

• Ensemble A 

Five original sets of experimental data as compiled by Amaldi and Schubert at ^/s = 
23.5, 30.7, 44.7, 52.8, and 62.5 GeV 0g. 

• Ensemble B 

Five sets of ensemble A including at each energy the data at y/s = 27.5 GeV M. 
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We want to investigate the effects in the eikonal (transfer momentum space) coming from 
the fits to both ensembles. 

We make use of the following analytical parametrization for the scattering amplitude 
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where acj, j3j, j = 1, 2, n are real free parameters and p(s) is the experimental p value at 
each energy ||. With this amplitude, we analyze the differential cross section 

| = 7 r|F(g,,)| 2 , -t = q\ (2) 

Fits to experimental data from ensembles A and B were performed through the CERN- 



MINUIT routine, which gives also the errors Act,-, A/3j in the free parameters |T^]. In this 
procedure all the parameters were completely free, taking different values at each energy. 
With this, no assumption at all was made on the dependences of the free parameters with 
energy or momentum transfer, allowing a good statistical result: 1.0 < x 2 /Ndf < 2.1 for 
a number of between 116 and 226, in all cases ||14|| . Propagation of the errors to 



the differential cross section [|13| takes account of extrapolated curves, which, as stressed by 
Lombard Q , cannot be excluded on statistical grounds. The complet fit results for ensembles 
A and B are displayed in Fig. |l] |L4| (numerical values of the free parameters are available 
from the authors). 

The complex eikonal in the momentum transfer space, x{Qi s )i is connected with the 
scattering amplitude, Eq. (1), through the well-known formulas 

F(q, s) = ij bdbJ (q, b)T(b, s) = i(T(b, a)), (3) 
T(b,s) = l-e ix{b ' s \ (4) 
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X(q,s) = (x(b,s)), (5) 

where Jo is the Bessel function (azimuthal symmetry assumed), T(b, s) is the profile function 
and the angular brackets denote a symmetrical two-dimensional Fourier transform. 

Here we are only interested in the imaginary part of the eikonal, xii which from Eq. 
(impact parameter space) reads 

Xi(b,s) = In 1 =. (6) 

/r?(M + [i - rVMP 



The complex profile T R (b, s) +iTj(b, s) is analytically determined through Eqs. (1) and @, 
together with the corresponding errors by propagation. Since we found that, even taking 
account of the errors, 

r?(M) «1, (7) 



[i-r fl (MP 

the eikonal may be expressed by 

Xi{b,s) «ln- 1 (8) 

i - r R (b, s) 

From Eqs. (1) and (ff), this quantity may be evaluated in terms of the fit parameters and 
also the corresponding errors, A\i, by propagation. 

At last we should evaluate the eikonal in the momentum transfer space, Eq. (^|), and the 
errors Axi{q, s). However, due to the structure of our parametrization (1), the transform 
(H) cannot be analytically performed and so the errors cannot be estimated as done in the 
previous steps. We recall that this is a typical problem in current empirical analysis. 

In order to evaluate the erros Axi{q, s) we introduce the following method. Generically 
we can expand Eq. (|8|) in the form 

Xi{b,s)=T R {b,s)+D{b,s), (9) 

where D(b, s) corresponds to the remainder of the series. Performing the Fourier transform 
we obtain 
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X i{q,s) = F I {q,s) + D{q,s). (10) 

Since the amplitude Fj(q,s) and errors AFj(q,s) are directly given by the fits, our task 
concerns the evaluation of 

D(q,s) = (D{b,s)), (11) 

with the corresponding errors, AD(q, s), and this is the central point of the method. 
First, from Eqs. (§) and (§), the quantity D(b,s) may be evaluated 

D(b, s) = ln[ 1 ] - F R (b, s), (12) 
1 - T R (b, s) 

and also the errors, AD(b, s), through error propagation from r R (b, s). Next, making use of 
the CERN-MINUIT routine this set of points with errors, D(b, s) ± AD(b, s), was fitted by 
a sum of Gaussians 

D f (b,a) = j2Aie- B '». (13) 

With this parametrization, D(q,s) in Eq. ([□]) may be analytically evaluated and also the 
errors, AD(q,s), may be estimated through the propagation of the errors in Aj, Bj, as 
given by the routine. At last, Eq. fllPP leads to s) and the error propagation furnishes 
Ax/(g, s). 



This method was used by Furget in order to determine the eikonal Xi(<7) s ) [ 15| - A novel 
aspect of our approach is its use in the estimation of errors. 

Two typical results, at y/s = 30.7 and y/s = 52.8 GeV, are displayed in Figs. | and 
|3] for ensembles A and B, respectively. We plotted the quantity Xi(qi s ) multiplied by q 8 , 
q 2 = — t 0. In the case of ensemble A, experimental data at \/s = 30.7 GeV are available 
up to t max = 5.75 GeV 2 and from Fig. [2] we see that no zero can be inferred. However, from 
Fig. |, the data at large q 2 (ensemble B, with t max = 14.2 GeV), lead to statistical evidence 
(with the uncertainties) for the existence of zero. This effect of the data at large momentum 
transfer is also evident at ^/s = 52.8 GeV (Figs. Q and ^) since in ensemble A experimental 
data at this energy are availiable up to t max = 9.75 GeV 2 . 



We conclude that data at large momentum transfer allow the identification of eikonal 
zeros (change of sign). Quantitatively the position of the zero and its uncertainties may be 
extracted from the plots as shown in Figs. ^| and 0. The results from ensembles A and B 
are displayed in Fig. |j. Due to the uncertainties and, meanly, its asymmetry around some 
points, the dependence of the position of the zero with the energy is difficult to be inferred 
on statistical grounds. 

However, if we disregard the asymmetries and take the uncertainty at larger momentum 
transfer symmetrically around each point, we are lead to the following quantitative estimates. 
Parametrizing the results from ensemble B (Fig. f|(b)) by 

ql = A + B\ns (14) 

gives A = 10.66 ± 3.75 GeV 2 , B = -0.9917 ± 0.993 GeV 2 with x 2 /N DF = 1.27. Although 
this indicates some decreasing with the energy, the slope is nearly compatible with zero. 
Assuming zero slope we obtain q 2 . = 6.9 ± 0.7 GeV 2 with x 2 I^df = 2.07. 

By the other hand, taking account of all the uncertainties displayed in Fig. |](b), we 
could roughly estimate the position of the zero at 7.0 ± 2.0 GeV 2 in the interval 23.5 GeV 
< v/i < 62.5 GeV. 

Extensions of our analysis at ^Js = 13.8 and 19.4 GeV are in course and this will be 
crucial for extracting the dependence of the zero with the energy. Also, detailed discussions 
on physical interpretations, connections with nonperturbative QCD, phenomenological ap- 
proaches and comparisons with other analyses are being concluded. 
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FIGURES 

FIG. 1. Fits to pp differential cross section data through Eqs. (1) and (||) between 23.5 GeV 
(upper) and 62.5 GeV (down): (a) ensemble A; (b) ensemble B. Curves and data were multiplied 
by factors of 10 ±4 |§. 

FIG. 2. Extracted imaginary part of the eikonal multiplied by t 4 from ensemble A: (a) =30.7 
GeV; (b) ^ = 52.8 GeV. 

FIG. 3. Same as Fig. [2] from ensemble B. 

FIG. 4. Position of the zeros in the eikonal as function of the energy from ensemble A (a) and 
ensemble B (b). 
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